REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704-0188)  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

08-01-2008  Journal  Article 

1  3.  DATES  COVERED  f From  -  To) 

4.  TITLE  AND  SUBTITLE 

A  Real-Time  Nearshore  Wave  and  Current  Prediction  System 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

0602435N 

6.  AUTHOR(S) 

Richard  A.  Allard,  James  D.  Dykes,  Yuan-Huang  L.  Hsu,  James  M.  Kaihatu, 
Daniel  Conley 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

73-6646-05-5 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NRL/J  A/7320-05-5281 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

800  N.  Quincy  St. 

Arlington,  VA  22217-5660 

10.  SPONSOR/MONITOR  S  ACRONYM(S) 

ONR 

11.  SPONSOR/MONITOR’S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

A  nearshore  wave,  tide  and  current  prediction  system  was  demonstrated  during  the  MREA04  Trial  in  the  Portuguese  coastal  waters  near  Pinheiro  do 
Cruz  during  the  early  spring  of  2004.  Daily  forecasts  of  regional  scale  wave  and  tidal  information  and  nearshore  waves  and  currents  were  generated  in 
DIOPS  utilizing  a  suite  of  regional  and  nearshore  models  forced  with  data  from  meteorological  and  oceanographic  production  centers.  A  limited  beach 
experiment  was  conducted  with  three  Nortek  current  meters  deployed  in  the  surf  zone  and  a  video  imagery  system  to  generate  10-min  time  exposures 
used  to  identify  the  locations  of  wave  breaking.  In  this  study,  Delft3D,  a  coastal  hydrodynamic  modeling  system,  capable  of  simulating  hydrodynamic 
processes  due  to  waves,  tides,  rivers,  winds  and  coastal  currents,  is  used  to  predict  the  nearshore  wave  and  longshore  current  near  Pinheiro  de  Cruz. 
The  nearshore  bathymetry  used  in  this  study  was  based  on  LIDAR  data  collected  in  February  2000.  Delft3D  shows  mixed  results  when  compared  with 
the  measured  wave  height  and  nearshore  currents.  Improved  Delft3d  results  can  be  achieved  in  the  future  if  migrating  sand  bars  can  be  measured 
and  included  in  the  modeling. 


15.  SUBJECT  TERMS 

Atmospheric  forcing,  bathymetry,  longshore  currents,  surf  zone,  water  level,  wave  breaking 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

OF 

PAGES 

22 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Richard  A.  Allard 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Unclassified 

Unclassified 

Unclassified 

UL 

19b.  TELEPHONE  NUMBER  (Include  area  code) 
228-688-4894 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


ELSEVIER 


Available  online  at  www.sciencedirect.com 

•  "  *  * 

ScienceDirect 


Journal  of  Marine  Systems  69  (2008)  37-58 


JOURNA  L  C)  I 

MARINE 

SYSTEMS 


www.clscvicr.com/locatc/jmarsys 


A  real-time  nearshore  wave  and  current  prediction  system 

Richard  Allard  3'*,  James  Dykes3,  Y.L.  Hsua,  James  Kaihatu3,  Daniel  Conley'1 

Naval  Research  Laboratory,  Oceanography  Division,  Stennis  Space  Center.  MS.  United  States 
h  NA  TO  Undersea  Research  Centre,  La  Spezia,  Italy 

Received  6  July  2005;  accepted  2  May  2006 
Available  online  22  February  2007 


Abstract 

A  nearshore  wave,  tide  and  current  prediction  system  was  demonstrated  during  the  MREA04  Trial  in  the  Portuguese  coastal 
waters  near  Pinheiro  da  Cruz  during  the  early  spring  of  2004.  Daily  forecasts  of  regional  scale  wave  and  tidal  information  and 
nearshore  waves  and  currents  were  generated  in  DIOPS  utilizing  a  suite  of  regional  and  nearshore  models  forced  with  data  from 
meteorological  and  oceanographic  production  centers.  A  limited  beach  experiment  was  conducted  with  three  Nortek  current  meters 
deployed  in  the  surf  zone  and  a  video  imagery  system  to  generate  10-min  time  exposures  used  to  identify  the  locations  of  wave 
breaking.  In  this  study,  Delft3D,  a  coastal  hydrodynamic  modeling  system,  capable  of  simulating  hydrodynamic  processes  due  to 
waves,  tides,  n\ers,  winds  and  coastal  currents,  is  used  to  predict  the  nearshore  wave  and  longshore  currents  near  Pinheiro  da  Cruz. 
The  nearshore  bathymetry  used  in  this  study  was  based  on  LIDAR  data  collected  in  February  2000.  Delft3D  shows  mixed  results 
when  compared  with  the  measured  wave  height  and  nearshore  currents.  Improved  Delft3D  results  can  be  achieved  in  the  future  if 
migrating  sand  bars  can  be  measured  and  included  in  the  modeling. 

€/  2007  Elsevier  B.V.  All  rights  reserved. 
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I.  Introduction 

A  wave,  tide  and  surf  prediction  system  was 
implemented  during  the  2004  Maritime  Rapid  Environ¬ 
mental  Assessment  (MREA04)  Trial  for  the  Portuguese 
coastal  waters  southwest  of  Lisbon  in  late  March 
through  April  10,  2004.  Real-time  forecasts  of  nearshore 
wave,  tide  and  surf  conditions  for  beaches  near  Pinheiro 
da  Cruz,  Portugal,  were  generated  aboard  the  NRV 
Alliance.  Military  critical  operations  may  require  a 
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nested  modeling  approach  to  produce  reasonable 
estimates  (and  forecasts)  of  surf  conditions  for  specific 
locations  of  interest.  Inputs  into  these  models  include 
atmospheric  forcing  (winds),  bathymetry,  directional 
wave  spectra  and  water  levels.  How  good  are  the  inputs 
into  such  models?  Arc  the  boundary  conditions  more 
important  than  the  wind?  This  paper  describes  the 
D10PS/Dclft3D  modeling  approach,  environmental 
inputs  into  these  models  and  comparisons  with  in  situ 
data.  Although  some  of  the  input  data  were  limited,  we 
address  how  such  a  wave  forecasting  system  can 
perform  in  an  “operational”  setting.  The  numerical 
results  presented  in  this  paper  are  from  model  hindcasts 
performed  after  the  MREA04  Trial.  The  reasons  for  this 
arc  two- fold:  ( 1 )  sensitivity  studies  were  performed  w  ith 
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a  triple-nested  wave  model  forced  with  winds  from  three 
different  atmospheric  models,  and  (2)  computational 
resources  available  during  MREA04  did  not  permit 
real-time  high-resolution  nearshore  modeling  with  ade¬ 
quate  temporal  resolution  to  compare  against  in  situ 
observations. 

This  paper  is  divided  into  seven  sections.  Section  2 
provides  some  perspective  on  the  MREA04  beach  ex¬ 
periment  and  Rapid  Environmental  Assessment.  Section 
3  discusses  the  components  of  the  Distributed  Integrated 
Ocean  Prediction  System  (DIO PS).  In  Section  4  we 
describe  the  nearshore  modeling  system  Dclft3D. 
Section  5  describes  the  observation  and  measurement 
systems  used  in  this  study.  The  results  and  analysis  of 
these  studies  are  detailed  in  Section  6.  We  discuss  our 
findings  in  Section  7. 

2.  Background 

The  MREA04  Trial  provided  an  opportunity  to  test 
and  evaluate  a  beach  environmental  reconnaissance  and 
monitoring  system  for  amphibious  warfare  support.  The 
study  area  near  Pinheiro  da  Cruz  was  utilized  during  the 
NATO  Linked  Seas  2000  Exercise  (Allard  et  al.,  2000), 
an  area  used  routinely  for  naval  and  amphibious 
training.  The  limited  beach  experiment  was  conducted 
by  implementing  high-resolution  wave  and  surf  models 
to  depict  the  nearshore  environment  and  deploy  real- 
time  data  acquisition  systems  for  situational  analysis 
and  bathymetry  updates. 

The  beach  experiment  was  conducted  from  April  2 
to  11,  2004,  near  Pinheiro  da  Cruz.  Instrumentation 
included  video  cameras  and  a  high-resolution  digi¬ 
tal  camera,  three  very  shallow  water  Nortek  current 
meters  in  the  surf  zone  and  a  portable  beach  meteoro¬ 
logical  station.  The  systems  contained  data  logging 
and  real-time  data  transfer  systems  through  satellite 
communications. 

3.  Modeling  components 

3.1.  DIO  PS 

The  Distributed  Integrated  Ocean  Prediction  System 
(DIOPS)  is  a  wave,  tide  and  surf  prediction  system 
(Allard  ct  al.,  2005)  composed  of  the  SWAN  (acronym 
for  Simulating  Waves  Nearshore)  wave  model;  a  re¬ 
locatable  tide  model,  PCTides;  and  the  1-D  Navy 
Standard  Surf  Model  (NSSM).  SWAN  is  typically 
initialized  by  offshore  directional  wave  spectra  obtained 
from  an  operational  production  center.  During 
MREA04,  a  triple-nested  DIOPS  SWAN  forecast  was 


produced  daily,  with  the  innermost  100-m  resolution 
nest  providing  boundary'  conditions  to  a  very  high 
resolution  coupled  wave/hydrodynamic  modeling  sys¬ 
tem  called  Delft3D.  Details  on  Delft3I)  are  discussed  in 
Section  4. 1 . 

During  MREA04,  DIOPS  predictions  were  per¬ 
formed  on  a  400-MHz  Sun  Blade  100  UNIX  worksta¬ 
tion,  while  the  PC-based  Dclft3D  simulations  were 
executed  on  a  Dell  Latitude  C640  PC,  with  I  GB  of 
memory  and  a  2.1 -GHz  CPU.  In  the  summer  of  2004, 
DIOPS  was  ported  to  a  windows  based  environment, 
allowing  all  calculations  to  be  performed  in  the  same 
(and  faster)  PC  environment.  In  this  study,  while 
running  on  a  Windows-based  platform  with  a  3.1 -GHz 
CPU,  a  48-h  prediction  takes  approximately  3  min  for 
the  tide  model  PCTIDES,  54  min  for  SWAN  (host  and  2 
embedded  nests,  with  the  innermost  nests  utilizing 
PCTIDES  water  levels)  and  2  h  for  Del 03 D. 

The  Dynamic  Information  Architecture  System 
(DIAS),  developed  by  the  Argonne  National  Laboratory, 
is  a  flexible,  extensible,  object-oriented  framework  for 
developing  and  maintaining  complex  simulations.  DIAS 
is  the  DIOPS  software  object  that  allows  these  models  to 
share  needed  information.  The  object-oriented  DIOPS 
framework  decomposes  the  maritime  environment  from 
the  deep  ocean  to  the  shore  into  classes  of  software 
objects,  each  with  its  own  spatially  distributed  sets  of 
attributes,  and  with  dynamic  behaviors  that  are  imple¬ 
mented  by  the  appropriate  ocean  physics  models.  The 
models  (c.g.,  SWAN  and  PCTides)  communicate  only 
with  domain  objects  and  never  directly  with  each  other. 
This  makes  it  relatively  easy  to  add  or  swap  models 
without  recoding.  An  atmosphere  object  provides  the 
required  meteorological  forcing  for  all  the  DIOPS  model 
components.  Fig.  1  depicts  the  DIOPS  model  domain 
starting  in  deep  water  progressing  to  the  surf  zone.  The 
user  has  the  option  of  using  the  1-D  NSSM  or  Delft3D 
for  surf  prediction.  A  series  of  shell  scripts  processes 
output  from  SWAN  and  PCTides  and  generates  input 
files  for  Delft3D  calculations. 

A  DIOPS  scenario  consisting  of  SWAN  and  PCTides 
was  set-up  and  loaded  prior  to  the  MREA04  Trial.  Daily 
updates  of  wave  spectra  and  meteorological  forcing 
aboard  the  NRV  Alliance  were  required  to  generate  daily 
48-h  predictions  of  wave  and  water  levels. 

3.2.  Meteorological  models 

Atmospheric  model  predictions  available  to 
MREA04  Trial  participants  were  coordinated  through 
different  meteorological  or  defense  institutions  in 
Germany,  Portugal,  France  and  the  United  Stales.  The 
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l  ig.  1.  Schematic  diagram  showing  DIOPS  model  flow  from  deep  water  to  the  surf  /one. 


following  section  briefly  describes  the  three  models 
used  in  this  study. 

The  U.S.  Navy's  Coupled  Ocean/Atmosphere  Me- 
soscalc  Prediction  System  (COAMPS&1 )  has  been 
developed  by  the  Naval  Research  Laboratory  (Hodur, 
1997)  and  is  run  operationally  at  the  Fleet  Numerical 
Meteorology  and  Oceanography  Center  (FNMOC)  for 
numerous  regions  around  the  globe.  CO  AM  PS  consists 
of  an  atmospheric  data  assimilation  system  comprising 
data  quality  control,  analysis,  initialization  and  non¬ 
hydrostatic  forecast  model  components.  The  horizontal 
resolution  of  the  COAMPS  Europe  application  grid 
used  in  this  study  was  0.2°.  Forecasts  were  available  at 
six-hourly  intervals  out  to  48  h.  COAMPS  Europe 
covers  the  entire  Mediterranean  originating  at  15°W. 
COAMPS  Europe  was  also  used  to  drive  the  WAVE- 
WATCH  111  wave  model  (Tolman,  2002)  that  provided 
deepwater  boundary  conditions  to  the  SWAN  model 


1  COAMPS  is  a  registered  trademark  of  the  Na\al  Research 
Laboratory. 


described  in  the  next  section.  WAVE  WATCH  111  is  a 
third  generation  wave  model  developed  at  NOAA/ 
NWS/NCEP. 

For  the  present  study,  a  particular  set-up  of  the 
METEO  FRANCE  ALAD1N  (Aire  Limitee  Adaptation 
dynamique  Dcveloppcmcnt  InterNational)  model 
(http://www.cnmi.mctco.fr/aladin)  was  used  and  pro¬ 
vided  by  the  French  Service  Hydrographique  et  Ocean- 
ographique  de  la  Marine.  The  model  domain  extended 
from  1 1  to  8°W  and  from  36  to  40°N  with  a  spatial 
resolution  of  0.1°  with  forecast  increments  of  3  h 
available  to  54  h  twice  daily. 

The  Deutscher  Wcttcrdicnst  Lokal-Modcll  (DWD- 
LM)  (http://www.dwd.dc.en)  covered  the  area  from  34 
to  42°N  and  15  to  6°W.  DWD-LM  received  boundary 
conditions  from  the  DWD  global  model.  The  non¬ 
hydrostatic  DWD-LM  performs  a  continuous  data 
assimilation  based  on  an  observational  nudging  scheme 
(Dorns  and  Schattler,  1999).  The  LM  is  run  operation¬ 
ally  with  a  7-km  horizontal  mesh  and  35  vertical  layers. 
The  data  provided  during  MREA04  were  interpolated  to 
a  0.05°  horizontal  resolution  grid.  DWD-LM  contained 
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f  ig.  2.  Domain  for  triple-nest  SWAN  wave  model  computations 
performed  during  MR1:A04.  Black  circles  indicate  locations  of 
WAV H WATCH  III  directional  wave  spectra  applied  as  boundary 
conditions  to  SWAN  host  grid.  Solid  black  rectangles  indicate  SWAN 
nests.  Black  star  denotes  locations  of  wave  buoy  near  Sines. 

the  highest  resolution  of  the  three  models  examined  in 
this  study. 

3.3.  SHAN 

SWAN  (http://fluidmechanics.tudelfl.nl/swan)  is  a 
third-generation,  phased-averaged  wave  model  (Booij 
et  al.,  1999;  Ris  ct  al.,  1999)  applicable  at  any  scale,  but 
most  efficient  when  predicting  wave  conditions  for  small 
scales.  It  is  capable  of  modeling  coastal  regions  with 
shallow  water,  barrier  islands,  tidal  flats,  local  winds  and 
ambient  currents.  SWAN  is  based  on  the  spectral  action 
balance  equation,  treated  in  discrete  form.  Short-crested, 
random  wave  fields  propagating  simultaneously  from 
widely  varying  directions  can  be  accommodated.  The 
SWAN  model  accounts  for  shoaling,  refraction,  wave 
generation  due  to  w  ind,  energy  dissipation  due  to  white¬ 
capping,  bottom  friction,  depth- induced  breaking  and 

Tabic  1 

Information  about  SWAN  model  configuration  used  during  MRHA04 


non-linear  wave  wave  interactions  (quadruplets  and 
triads).  Although  the  SWAN  version  used  in  these 
simulations  (40.11)  does  not  account  for  diffraction, 
version  40.41  has  included  diffraction  in  a  phase- 
decoupled  approximation. 

A  triple-nested  SWAN  wave  forecast  capability  was 
established  for  the  Portuguese  coastal  waters  for  MREA04; 
the  nests  are  shown  in  Fig.  2.  The  host  grid  extended  from 
36  to  40°N  and  1 1  to  8°W,  with  a  grid  resolution  of  0.05°. 
All  three  nests  contained  25  frequency  bins  ranging  from 
0.05  to  1 .0  Hz,  and  directional  bins  ranging  from  0  to  360° 
at  10°  intervals.  Boundary  conditions  of  directional  wave 
spectra  were  specified  on  the  outer  boundary  of  the  SWAN 
host  grid.  Embedded  nests  with  resolutions  of  0.01°  and 
0.001  °  were  contained  within  the  host  grid  with  boundary 
conditions  provided  via  SWAN  nesting  files.  Tabic  1 
provides  additional  information  about  grid  resolutions 
and  settings  used  for  these  simulations.  The  inner  nest  of 
SWAN  utilized  water  levels  from  PCTides  (sec  Section 
3.4)  to  adjust  the  water  depths  hourly. 

The  NATO  Undersea  Research  Centre  provided 
bathymetry  for  the  outer  two  SWAN  grids  at  resolutions 
of  6  and  12  s.  The  innermost  SWAN  nest  covered  an  1 8- 
km  section  of  the  beach  area  centered  near  Pinheiro  da 
Cruz.  The  bathymetry  for  this  1 00-m  resolution  nest  was 
based  in  part  on  the  SHOALS  (acronym  for  Scanning 
Hydrographic  Operational  Airborne  LIDAR  Survey) 
hydrographic  data  collected  during  the  Linked  Seas 
2000  NATO  Exercise.  SHOALS  w'as  deployed  on  a 
Twin  Otter  aircraft  in  April  2000  over  a  two-day  period 
to  collect  bathymetry  for  a  16-km  section  of  coastline 
ranging  from  depths  of  14  m  above  the  low  water  line  to 
depths  of  22  m  (Lillycrop  et  al.,  2000). 

The  southwest  coast  of  Portugal  is  predominantly 
affected  by  waves  propagating  from  the  northwest  and 
west,  with  occasional  motion  from  the  southwest.  Buoy 
statistics  from  the  Sines  buoy  (http://www.hidrografico. 
pt/wwwbd/Boias/BoiasUltimoRegisto.asp)  used  in  this 
study  indicate  a  predominant  wave  direction  from  the 
northwest  60%  of  the  time  for  the  one-year  period 
ending  April  1, 2005,  based  on  more  than  30,000  hourly 
observations.  The  Pinheiro  da  Cruz  beach  is  a  sandy- 
cliff  stretch  (Pires-Silva  et  al.,  2002)  nearly  centered  in  a 
soft  cell  coast  (approximately  60  km  in  length)  bounded 
by  the  Sado  Estuary  to  the  north  and  Sines  I  \ arbor  to  the 


No.  of  A'  No.  of  Y  LONGW  LONGK  LATS  LATN  RHS  No.  of  direction  No.  of  frequency  BC  Timestep 

Grid  A  61  SI  349  352  36  40  0.05  36  24  W'W'3  15  min  NON 

Grid  B  52  59  350.82  351.33  38  38.58  0.01  36  24  Grid  A  nest  15  min  NON 

Grid  C  81  181  351.16  351.24  38.15  38.33  0.001  36  24  Grid  B  nest  3-hourly  STAT 
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Fig.  3.  Bathymetry  near  Pinheiro  da  Cruz,  isobaths  displayed  at  50-m 
intervals.  Circle  near  coast  shows  beach  location;  triangle  denotes 
location  of  Meteo  buoy.  Black  outer  rectangle  denotes  SWAN  inner 
nest  (GRID  C).  Inner  rectangle  denotes  DeIft3D  WAVE  grid. 


south.  The  bathymetry  (Fig.  3)  is  fairly  regular  with 
contours  paralleling  the  coast  with  the  exception  of  the 
Setubal  C  anyon  which  terminates  approximately  10  km 
northwest  of  Pinheiro  da  Cruz,  the  location  of  which  is 
indicated  by  the  black  circle  on  the  coastline  in  Fig.  3. 
During  events  with  strong  northerly  to  north- north¬ 
westerly  winds,  the  region  is  sheltered  from  higher  wave 
heights  due  to  the  short  fetch. 

The  SWAN  host  was  forced  with  boundary  conditions 
from  the  27-km  Europe  wave  model  WAVEWATCH  111 
provided  by  the  FNMOC.  Directional  wave  spectra  were 
provided  via  a  METCAST  server  for  the  locations 
specified  in  Fig.  2  throughout  the  MREA04  Trial. 
WAVEWATCH  111  data  were  not  available  during  the  00 
GMT  (Greenwich  Mean  Time)  forecast  cycles  on  March 
30  and  April  X,  2004.  Missing  fields  were  replaced  with 
forecast  data  from  the  previous  watch  cycle.  WAVE- 
WATCH  III  is  forced  with  winds  from  the  27-km 
COAMPS  Europe  model  runs. 

3.4.  PCTides 


components:  ( 1 )  a  grid  generator  that  utilizes  the  Naval 
Research  Laboratory's  DBDB2  global  bathymetry'  data¬ 
base  and  additional  higher  resolution  regional  databases; 
(2)  a  global  tidal  model,  the  Finite  Element  Solutions  99 
(FES99)  which  is  used  to  provide  tidal  conditions  on 
open  boundaries;  and  (3)  a  two-dimensional  barotropic 
ocean  model  that  produces  predictions  of  water  level 

Detft3d  Domains 
%  10®  UTM  Coordinate  System 


Easting  Coordinates  (m)  -»  „ 


PC' Tides  is  a  globally  re-locatable  tide/surge  forecast 
system  (Prcller  et  al.,  2002)  consisting  of  three  primary 


Fig.  4.  Domains  for  WAVE  denoted  by  large  rectangle  and  FLOW 
denoted  by  visible  bathymetry,  whose  depth  values  in  meters  are 
indicated  by  the  color  bar. 
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and  depth-averaged  currents.  Options  include  the  ability 
to  assimilate  International  Hydrographic  Office  (IHO) 
tidal  station  data,  wetting  and  drying  and  a  surge  capa¬ 
bility  that  utilizes  wind  forcing  and  sea  level  pressure 
from  operational  meteorological  models. 

During  MREA04,  PCTides  48-h  forecasts  of  water 
levels  were  provided  daily  for  a  grid  encompassing  the 
area  from  34  to  40°N  and  1 2  to  5° W  with  a  resolution  of 
0.05°  (PCTides  can  be  run  on  cither  a  Lambert  Con¬ 
formal  or  spherical  grid).  Tidal  constituents  for  the  M2, 
S2,  N2,  K.2,  Kl,  Ol,  PL  Ql  and  2N2  tides  were  ex¬ 
tracted  from  the  Finite  Element  Solutions  99  (Lefevre 
et  al.,  2002).  DWD-LM  winds  and  sea  level  pressure 
were  used  in  the  daily  DIOPS  run  stream.  The  output 
was  fed  to  the  inner  nest  of  SWAN  to  adjust  bathymetry 
at  hourly  intervals.  A  total  of  28  IHO  stations  were 
assimilated  into  the  PCTides  model  solutions. 

4.  Near  shore  modeling  system 

4.1.  Modeling  description 

The  Del  03  D  system,  developed  by  Delft  Hydraulics 
(http://www.wldclft.nl),  is  a  complete  coastal  hydrody¬ 


namic  modeling  system,  capable  of  simulating  hydrody¬ 
namic  processes  due  to  waves,  tides,  rivers,  winds  and 
coastal  currents;  the  present  application  of  the  model  is 
focused  on  nearshore  hydrodynamics  forced  by  breaking 
waves.  The  model  can  be  run  in  rectangular  (equidistant 
or  stretched)  or  curvilinear  coordinates;  all  necessary  grid 
generation  software  for  creating  curvilinear  grids  is 
included  with  the  Dclft3D  package. 

The  Delft3D  system  uses  two  modules  for  simulating 
nearshore  wave-induced  hydrodynamic  processes.  The 
WAVE  module  uses  the  SWAN  model  for  propagation 
and  generation  of  waves.  Hydrodynamics  are  simulated 
with  the  FLOW  module  (WL  Delft  Hydraulics,  2001), 
which  uses  the  three-dimensional  hydrostatic  shallow 
water  Navier-Stokes  equations  with  the  discretization 
scheme  of  Stclling  and  Van  Kcster  (1094).  Forcing  for 
wave-induced  flow  is  derived  from  radiation  stresses 
with  both  wave  and  roller  contributions  (Renicrs  ct  al., 
2002;  Roelvink,  2003).  The  model  can  be  run  with  both 
one-way  forcing  or  with  feedback  between  the  two 
modules.  Although  a  highly  flexible  tool  for  various 
applications  small  and  large,  this  component  of  the 
nearshore  modeling  system  was  tailored  specifically  for 
a  domain  that  would  extend  from  the  shoreline  to  about 
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I  ig.  5.  Location  of  surf  sensors  at  Pinheiro  da  Cruz.  Sensors  arc  overlain  on  rectified  time  exposed  image  of  the  surf  zone  created  at  low  tide  on  April  7. 
The  light  areas  in  the  image  represent  regions  of persistent  wave  breaking  and  indicate  the  bar  locations.  Model  results  arc  compared  primarily  w  ith  the 
measurements  collected  from  the  Nortek  velocimeler  N6. 
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I  km  seaward  (Dykes  ct  al.,  2003).  Morris  (2001 )  first 
implemented  and  evaluated  Delft3D  with  this  specific 
application  in  mind. 

Originally  the  NSSM  provided  surf  conditions  for  a 
one-dimensional  beach  profile  across  the  shore.  Given 
accurate  inputs  from  the  host  SWAN  model,  predictions 
from  this  model  were  sufficiently  accurate  provided 
there  was  minimal  depth  variability  along  the  shore. 
However,  when  sufficient  variability  in  water  depths 


exists  with  bars  and  bar  breaks,  the  one-dimensional 
predictions  cannot  account  for  currents  due  to  the 
resulting  variable  longshore  pressure  gradient  forces. 
Additionally,  the  NSSM  can  be  considered  a  highly 
constrained  model,  in  the  sense  that  wave  breaking  and 
dissipation  drives  only  the  longshore  current,  with  no 
ability  to  generate  cross-shore  flows  like  rip  currents. 

In  this  exercise,  we  used  the  Dclft3D  WAVE  and  two- 
dimensional  FLOW  modules.  WAVE  uses  the  SWAN 


Month/Day  2004 

Fig.  6.  (A)  Sea  level  pressure  (mb)  from  COAMPS,  ALADIN,  DWD-LM  and  Mctco  buoy  for  the  period  April  4  11;  (B)  same  as  a.  but  for  wind 
direction.  Buoy  data  are  missing  from  April  4  9  ( IS  GMT);  (C)  same  as  a,  but  for  wind  magnitude  (m/s). 
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Fig.  6  (  continued ). 


wave  model  in  stationary  mode.  A  master  module  coor¬ 
dinates  the  two  components  which  are  coupled  together 
to  simulate  wave  current  interaction.  A  communication 
file  acts  as  the  interacting  transfer  staging  point  for 
exchanging  data  between  the  two  models. 

4.2.  Domain  and  bathymetry 

Using  RGFGR1D,  a  grid-building  tool  of  Delft3D, 
domains  for  the  WAVE  and  FLOW  modules  were 
established  based  on  the  availability  of  the  water  depth 
data  and  the  location  of  the  Pinheiro  da  Cruz  beach.  Due 
to  their  geographic  orientation,  both  computational  grids 
could  be  laid  conveniently  into  the  Universal  Transverse 
Mercator  (UTM)  coordinate  system.  The  resolution  of 
the  WAVE  domain  was  about  20  m  in  the  ^-direction 
and  40  m  in  the  F-direction,  while  that  of  FLOW  varied 
from  10  to  30  m.  On  each  side  of  the  WAVE  domain,  the 
lateral  boundaries  were  placed  about  500  m  beyond  the 
FLOW  lateral  boundaries  away  from  the  center  to  keep 
the  wave  forcing  from  introducing  spurious  forcing 
gradients  at  those  boundaries.  Thus  the  length  of  the 
FLOW  domain  was  about  2.5  km  along  the  shore  while 
the  WAVE  domain  turned  out  to  be  about  3.5  km.  The 
width  from  the  shore  of  both  domains  was  about  1  km. 
Fig.  4  depicts  the  relationship  between  the  two  grids. 
The  outer  box  is  the  WAVE  grid  and  the  FLOW  grid  is 
represented  by  the  coverage  of  bathymetry. 

Water  depths  taken  from  the  SHOALS  L1DAR  data 
collected  during  the  Linked  Seas  2000  exercise  were 


used  as  the  initial  depths  for  Delf\3D.  Using  the  Del 03 D 
tool  called  QIJ1CK1N  for  applying  bathymetry  to  the 
grid,  the  sampled  water  depths  were  interpolated  to  the 
WAVE  and  FLOW  grids.  The  sample  depths  were  dense 
enough  that  a  simple  grid  cell  averaging  routine  in 
QU1CKIN  was  sufficient  to  populate  all  the  grid  points 
with  some  minor  fill  in.  Fig.  4  depicts  the  water  depths 
used  for  the  FLOW  domain.  Within  about  500  m  of  both 
lateral  boundaries  of  the  FLOW  domain  are  artificially 
extended  values  to  allow  for  a  smooth  transition  from 
the  boundaries.  The  purpose  of  this  smooth  transition  is 
to  prevent  the  lateral  boundaries  from  generating 
artificial  secondary  flows  that  may  enter  the  domain 
(Roelvink  and  Walstra,  2004).  Values  not  shown  for  the 
WAVE  domain  are  also  extended  for  the  same  reason. 

Initially,  all  the  grid  points  in  the  FLOW  and  WAVE 
grids  were  equidistant  and  perfectly  rectilinear.  To 
enhance  computational  performance,  the  grid  points  for 
the  FLOW  domain  were  later  thinned  out  at  points  more 
distant  from  the  area  of  interest. 

In  addition  to  accounting  for  a  two-dimensional 
domain  of  variable  water  depths,  Delft3D  can  incorpo¬ 
rate  other  two-dimensional  inhomogeneous  boundary 
conditions  such  as  winds,  water  levels  and  wave  spectra. 
In  this  case,  the  domain  was  deemed  small  enough  to 
assume  a  uniform  wind  and  water  level.  However,  at  the 
boundaries  there  were  spatially  varying  wave  conditions 
which  were  delivered  in  the  form  of  spectra  data  output 
from  the  SWAN  innermost  nest  (Grid  C)  in  DIOPS, 
giving  six  points  for  each  lateral  side  and  ten  points  on 
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the  seaward  side.  This  level  of  detail  may  be  important 
particularly  when  bars  arc  not  present  as  the  variable 
wave  energy  may  still  produce  rips. 

All  the  details  associated  with  the  domain  grids, 
water  depths,  the  numerical  and  physical  parameters 
and  outputs  were  set-up  ahead  of  time  using  the 


Delft3D  graphical  user  interface.  Except  for  some  small 
adjustments  in  the  numerical  behavior,  no  changes  were 
needed  during  the  exercise.  In  an  operational  context, 
updated  bathymetry  can  be  incorporated  into  the 
models  without  interfering  with  the  normal  course  of 
operations. 
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t  ig  7.  10-m  winds  from  (a)  COAMPS,  (b)  ALADIN  and  (c)  DWD-LM  on  April  6  at  0,  6.  12  and  IS  CiMT.  Data  arc  plotted  al  0.2°  intervals 
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Table  2 

Meteorological  model  specifics 


LONCiW 

LONGK 

LATS 

LATN 

RES 

Output 

frequency 

(h) 

Model 

ALADIN 

349 

352 

36 

40 

0.1° 

6 

DWD-LM 

345 

354 

34 

42 

0.05° 

3 

COAMPS 

345 

40 

29 

59 

0.2° 

3 

4.3.  MREA04  nearshore  forecast  operations 

Each  day  DIOPS  produced  a  48-h  forecast  starting 
at  00  GMT,  resulting  in  a  set  of  files  with  spectra, 
winds  and  water  levels  retained  for  Delft3D.  Utilities 
were  developed  to  extract  the  needed  information  from 
those  files  and  to  incorporate  that  information  into 
the  Dclft3D  configuration  files,  all  in  an  automated 
fashion. 

Initially,  due  to  some  time  dependence  limitations  in 
the  WAVE  module,  the  run-time  strategy  was  to  run 
individual  simulations  of  WAVE  and  FLOW  for  certain 
points  in  time.  The  SWAN  model  in  the  WAVE  module 
provided  forcing  for  FLOW  which  would  run  for 
60  min,  the  time  period  considered  sufficient  to  reach 
equilibrium  holding  all  other  conditions  such  as  water 
level  and  wind  constant.  This  configuration  in  effect 
resulted  in  stationary  conditions,  i.c.,  to  a  point  of 
equilibrium.  In  this  way,  we  chose  to  run  Dclft3D  for 
every  3  h  throughout  the  forecast  period.  Coupling  to  the 
wave  current  interaction  w'as  one-way  (waves  forcing 
currents). 

After  the  exercise,  a  time-dependent  capability  was 
incorporated  into  the  WAVE  module  software  allowing 
for  model  suite  execution  to  be  performed  in  a  con¬ 
tinuous  fashion  after  a  one-time  set-up.  For  a  hindcast  of 
the  entire  period,  FLOW  ran  continuously  interrupted 
only  by  the  hourly  mil  of  the  WAVE  module.  This  mode 
allowed  for  the  input  of  time  series  of  water  levels  and 
wind  inputs.  The  WAVE  module  can  interpolate  from 
the  spectra  input  and  can  keep  its  continuity  using  restart 
files.  Two-way  coupling  to  simulate  wave-current 
interaction  in  a  real-time  situation  is  now  feasible. 

5.  Observations  and  measurements 

5.1.  Sutf  zone  wave  and  current  measurements 

Measurements  of  waves  and  currents  in  the  nearshore 
were  collected  by  deploying  three  different  Nortek 
Vector  acoustic  vclocimcters  in  the  surf' zone  (Fig.  5). 


These  sensors  arc  capable  of  simultaneously  measuring 
the  fluid  pressure  and  three  orthogonal  components  of 
velocity  at  a  single  location.  The  sensors  were  equally 
spaced  in  the  longshore  spanning  a  total  length  of 
100  m.  At  the  time  of  installation,  all  sensors  were 
placed  at  low  tide  in  1.5  m  w'ater  depth  in  the  shoreward 
most  trough  of  the  study  site.  The  two  southern  most 
sensors  (N9  and  N6)  were  connected  to  a  shore-based 
laboratory  by  means  of  a  data  cable  through  which 
power  was  supplied  to  the  sensors  so  that  two-way 
communication  permitted  remote  control  of  the  data 
collection  procedure  and  real-time  data  telemetry.  The 
raw'  data  were  collected  at  8  Hz  and  stored  for  later 
processing.  Archived  data  were  cleaned  as  described  by 
Elgar  et  al.  (2001)  and  reduced  to  2  Hz.  The  north¬ 
ernmost  sensor  (N4)  was  developed  as  part  of  an 
autonomous  beach  monitoring  system  and  data  proces¬ 
sing  occurred  in  the  accompanying  underwater  package 
after  which  the  analyzed  parameters  w  ere  transmitted  to 
a  home  laboratory  through  Iridium  satellite  modem.  For 
this  sensor,  analysis  was  based  on  collection  of  12  min 
of  8  Hz  data  every  2  h. 

5.2.  Imagery  at  the  beach 

An  image  monitoring  system  for  the  surf  zone  was 
also  deployed  as  part  of  the  MREA04  experiment.  This 
system  contained  both  a  digital  camera  and  color 
digital  video  camera.  The  digital  camera  w'as  used  to 
collect  high-resolution  (3.7  M  pixel)  snapshots  of  the 
site  and  the  video  camera  was  used  to  collect  high- 
frequency  (^2  Hz)  images  for  a  sustained  period  of 
time  from  w'hich  a  single  time  exposed  image  (Lippmann 
and  Holman,  1989)  was  formed.  Both  cameras  were 
mounted  on  a  common  pan-tilt  device  and  images  span¬ 
ning  the  entire  beach  were  collected  by  progressively 
repositioning  the  field  of  view  of  the  cameras  from  south 
to  north.  The  cameras  were  mounted  on  the  bluff  behind 
the  beach  at  an  altitude  of  22  m  and  a  distance  of 
approximately  95  m  behind  the  shoreline  (Fig.  5).  The 
monitoring  system  was  programmed  to  collect  a 
complete  series  of  four  images  from  each  camera  three 
times  a  day  (morning,  noon  and  evening).  The  time- 
averaged  images  were  processed  onboard  the  system 
central  control  unit  and  the  final  images  were  transmit¬ 
ted  remotely  using  FTP  protocol  over  a  V-Sat  Internet 
connection.  Using  techniques  (Holland  et  al.,  1 997)  which 
were  provided  by  the  ARGUS  program  (Holman  et  al.. 
1993;  http://cil-www.coas.oregonstate.cdu:8080/),  the 
images  were  subsequently  rectified  and  merged  to  form 
a  single  view7  of  the  study  site  (for  examples,  see  Figs.  5 
and  13). 
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6.  Results  and  analysis 

6.1.  Meteorological  models  versus  Meteo  buoy 

To  gain  a  better  understanding  of  the  deepwater  wave 
predictions  (next  section)  that  impact  the  nearshore 
models,  we  assess  the  performance  of  the  three  mete¬ 
orological  models.  During  MREA04,  an  Aanderaa 
Instruments  Coastal  Monitoring  Buoy  3280  (hereafter 


referred  to  as  the  Meteo  buoy)  was  deployed  west- 
southwest  of  Pinheiro  da  Cruz  at  a  120-m  water  depth. 
The  compact  data  buoy  is  designed  for  use  in  coastal 
waters,  ports  and  nearshore  platforms.  This  buoy  con¬ 
figuration  measured  wind  speed  (2.5-m  height),  wind 
direction,  sea  level  pressure,  wave  height,  air  temperature 
and  water  temperature.  Fig.  2  depicts  the  location  of  the 
Meteo  buoy.  The  buoy  data  were  not  assimilated  into  any 
of  the  meteorological  models.  The  processed  buoy  data 
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l;ig.  S.  (A)  Significant  wave  height  (meters)  during  the  period  March  31  April  11  (00  GMT)  at  the  Sines  buoy  location  versus  SWAN  forced  with 
ALADIN.  DWD-l.M  and  COAMPS  1 0-m  winds:  (B)  same  as  a.  but  for  mean  wave  direction  (°);  (C)  same  as  a,  but  for  peak  wave  period  (seconds). 
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were  available  at  approximate  one-hour  intervals.  In  the 
following  section,  three-hourly  buoy  data  arc  compared 
against  modeled  data  at  that  same  interval. 

Fig.  6A  shows  a  comparison  of  sea  level  pressure 
from  all  three  meteorological  models  versus  the  Meteo 
buoy  during  the  period  of  April  4-11, 2004.  T  here  is  very 
good  agreement  among  all  the  meteorological  models 
with  the  exception  of  higher  values  from  ALAD1N  on 
April  0.  Overall,  the  models  compare  well  against  the 
buoy,  although  the  buoy  consistently  reports  lower  values 
by  about  two  to  three  millibars  during  the  period. 

The  wind  directions  during  this  period  are  displayed 
in  Fig.  6B.  The  models  are  in  general  agreement 
throughout  the  period,  although  COAMPS  docs  not 
reflect  the  shift  from  offshore  to  onshore  flow  as  well  as 
the  ALADIN  and  DWD-LM  models.  This  could  be 
attributed  to  the  frequency  of  COAMPS  output  which 
was  available  at  six-hourly  interv  als  compared  to  three- 
hourly  for  the  other  models  in  addition  to  the  relatively 
coarse  27-km  horizontal  grid  resolution.  A  diurnal  cycle 


is  evident  throughout  the  period  as  heating  over  land 
generates  a  late  afternoon  sea  breeze.  The  pattern 
reverses  during  the  evening  hours.  Unfortunately,  the 
wind  direction  sensor  aboard  the  Meteo  buoy  was  not 
functioning  properly  until  the  afternoon  of  April  9. 
During  the  very  short  period  of  good  data,  there  is  very 
good  agreement  between  models  and  observation. 

The  Meteo  buoy  w  ind  speeds  were  adjusted  from  2.5 
to  10  m  using  the  simple  relationship  from  Hsu  et  al. 
(1994) 

U2  =  ul(:2/zi)'\  (1) 

where  u2  denotes  the  wind  speed  at  the  reference 
height  (2.5  m),  z2  and  it\  represents  the  wind  speed 
measured  at  height  Z\  (10  m).  The  exponent,  P,  is  set 
to  0.11  based  on  an  empirical  relationship  for  typical 
ocean  conditions. 

Fig.  6C  presents  a  comparison  of  modeled  10-m 
wind  speeds  versus  the  adjusted  10-m  wind  speed 


table  3 

SW'AN  wave  model  statistics  at  Sines,  Portugal,  forced  with  different  meteorological  models 

ALADIN  DWD-LM  COAMPS 


RMS 

Correlation 

Bias 

Scat 

RMS 

Correlation 

Bias 

Scat 

RMS 

Correlation 

Bias 

Scat 

Hs  (m) 

0.34 

0.94 

0.59 

0.1X6 

0.34 

0.96 

0.69 

0.189 

0.32 

0.95 

0.56 

0.177 

Direction  (°) 

11.7 

0.72 

127.2 

0.038 

10.08 

0.92 

-26.96 

0.033 

8.15 

0.94 

-45.21 

0.027 

TP  (s) 

1.35 

0.55 

6.28 

0.128 

1.05 

0.76 

3.66 

0.100 

1.09 

0.74 

3.98 

0.103 

R  Allan/  cl  al.  /Journal  of  Marine  Systems  69  (2008)  37  58 


49 


Fig.  9.  Ratio  of  significant  wave  heights  for  simulations  performed  with  three  different  wind  models  without  spectra  applied  on  the  model  boundary 
div  ided  by  the  observed  wave  height  at  the  Sines  buoy.  Low  percentages  indicate  that  swell  energy  dominates  while  high  values  suggest  that  local 
wind  generation  is  the  primary'  mechanism  for  wave  growth. 


observed  at  the  Meteo  buoy.  Overall,  the  observed  wind 
speeds  are  lower  than  modeled  by  more  than  2.0  m/s 
during  this  period.  We  hypothesize  that  an  over- 
prediction  in  wind  speed  by  the  numerical  models  may 
be  partially  responsible  for  the  higher  wave  heights  (see 
next  section)  indicated  by  SWAN  at  this  location. 

A  comparison  of  10-m  winds  for  the  region  bounded 
from  36  to  40°N  and  1 1  to  8°W  for  April  6  is  shown  in 
Fig.  7.  Data  are  plotted  on  a  0.2°  grid  to  account  for  the 
differing  model  resolutions  (DWD-LM  had  the  highest 
resolution  at  0.05°,  while  COAMPS  had  the  lowest 
resolution  at  0.2°).  While  all  three  models  indicate  a 
switch  from  north  to  north  northeasterly  How  to  north 
northwesterly  flow  by  1 8  GMT  near  the  Meteo  buoy,  the 
magnitude  is  weakest  from  CO  AM  PS. 

6.2.  SWAN  versus  SINES  buoy 

During  MREA04,  SWAN  predictions  were  made 
daily  aboard  the  NRV  Alliance  for  the  period  April  2- 
10,  2004.  Wind  forcing  consisted  of  the  DWD-LM 
meteorological  model,  covering  the  same  domain  as  the 
SWAN  host  grid,  with  a  resolution  of  0.05°.  In  this 
study,  we  performed  three  hindcasts  using  the  follow  ing 
meteorological  atmospheric  models:  (I)  DWD-LM,  (2) 
ALAD1N  and  (3)  COAMPS-Europe.  Table  2  provides 
information  about  these  models  including  horizontal 
and  temporal  resolution  of  forecast  fields. 


A  directional  deepwater  “WAVECM  buoy  moored  at  a 
depth  of  97  m  near  Sines,  Portugal  (denoted  as  a  star  in 
Fig.  2)  w'as  used  to  evaluate  the  performance  of  SWAN 
using  the  three  meteorological  models  described  earlier. 
Fig.  8A-C  show  a  comparison  of  significant  wave 
height,  mean  wave  direction  and  peak  wave  period 
during  the  period  March  31 -April  1 1,  2004.  Buoy  data 
were  available  at  four-hourly  intervals.  A  major  storm 
event  occurred  during  April  1-3  with  wave  heights  at 
the  buoy  exceeding  4.0  m  on  April  3.  An  intense  storm 
located  near  46°N,  13°W  on  April  1  produced  swells 
that  reached  this  area  a  few  days  later. 

Table  3  summarizes  the  comparisons  between  SWAN 
(host  grid)  at  the  Sines  buoy  forced  w'ith  the  three 
meteorological  models.  However,  since  the  directional 
wave  spectra  provided  by  WAVE  WATCH  Ill  is  forced 
with  COAMPS,  it  was  not  possible  to  determine  if  those 
other  models  (e.g.,  ALAD1N)  would  have  produced 
significantly  different  results.  Significant  wave  height 
RMS  error  statistics  show  that  all  three  models  per¬ 
formed  similarly  with  a  slight  edge  given  to  COAMPS 
whose  RMSE  ranged  between  32  and  34  cm  with 
correlation  coefficients  falling  between  0.94  and  0.96. 
Scatter  indices  for  all  three  test  cases  were  consistently 
in  the  1 8-  to  1 9-cm  range.  The  scatter  index  is  defined  as 
RMSE  divided  by  mean  of  the  observations. 

Fig.  8B  shows  that  the  hindcasts  performed  w  ith 
ALADIN  missed  the  turning  of  the  waves  from 
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northwesterly  to  southwesterly  on  April  1.  The  RMS 
wave  direction  errors  were  near  10°,  with  COAMPS 
displaying  the  lowest  error  (8.1°)  and  highest  correlation 
(/*= 0.94).  All  three  models  under-predicted  the  14-s 
peak  periods  shown  in  Fig.  8C  on  April  3  and  over- 
predicted  the  peak  period  on  April  10.  The  DWD-LM 


and  COAMPS  showed  RMS  errors  at  1.1  s,  with  the 
highest  correlation  (/*=  0.76)  with  the  DWD-LM. 
Overall,  all  three  meteorological  models  with  different 
resolutions  and  physics  produced  similar  results  when 
forcing  SWAN.  Scatter  indexes  for  peak  period  and 
wave  direction  were  small  for  all  three  test  cases. 


/\  Water  Levels  at  Sines 

PCTides  vs  Buoy  Scatterplot  1-11  April  2004 
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Time  Series:  PCTides  Model  vs  N6  Measurements 
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Fig.  10.  (A)  PCTides  water  levels  (meters)  versus  data  collected  by  FS  Borda;  (B)  PCTides  water  levels  (meters)  versus  data  collected  from  N6  near 
Pinheiro  da  Cm/;  (C)  Scatter  plot  of  PCTides  water  levels  versus  N6  data. 
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C  Scatterplot  PCTides  vs  N6  Measurements 


Measured  Water  Level  (m) 


Fig.  10  ((  ontinued ). 


The  hindcast  experiments  described  earlier  were  all 
driven  with  WAVEWATCH  111  wave  spectra  forced  with 
COAMPS-Europe  wind  fields.  We  did  not  have  wind 
forcing  available  for  the  Atlantic  Basin  to  examine  the 
performance  of  WAVEWATCH  III  with  the  other 
meteorological  models  or  their  operational  counterparts 
(c.g.,  DWD  global  model).  To  investigate  the  effect  of 
offshore  boundary  conditions  on  these  hindcasts,  three 
additional  experiments  were  performed  in  which 
boundary  forcing  was  not  applied  to  the  SWAN  host 
grid.  Each  experiment  was  forced  with  wind  only  from 
the  three  meteorological  models  described  earlier.  Fig.  9 
depicts  the  percentage  of  the  significant  wave  height 
(with  no  boundary'  forcing)  divided  by  the  observed 
wave  height  observed  at  the  Sines  buoy  location.  A 
value  of  100%  indicates  wind-generated  waves  are 
generated  without  regard  to  swell  generated  outside  of 
the  model  domain.  Lower  values  (e.g.,  10-50%)  imply 
that  the  propagation  of  swell  events  is  not  represented  by 
the  exclusion  of  boundary  forcing.  All  thr  ee  wind  models 
produced  similar  SWAN  wave  heights  during  the  period 
of  April  3-8  with  differences  between  all  three  tests  at 
less  than  20%.  After  that  period,  CO  AM  PS  and  DWD- 
LM  winds  produced  wave  heights  10-30%  higher  than 
observed  on  April  9,  while  ALADIN  wave  heights  show 
much  better  agreement  with  observ  ations.  To  summarize: 
when  energy  (e.g.,  from  WAVEWATCH  III)  from  outside 
the  SWAN  host  domain  was  dominant,  wave  heights  for 
all  three  hindcasts  were  similar.  However,  when  wind 
forcing  was  dominant  (c.g.,  on  April  2  and  10),  sig¬ 


nificant  differences  among  the  hindcasts  performed  with 
the  three  meteorological  models  was  evident.  Further 
review  into  this  topic  could  address  how  to  improve 
operational  w  ave  forecasts  where  local  wind  generation  is 
dominant. 

6.3.  SWAN  versus  Meteo  buoy 

The  Meteo  buoy  was  damaged  during  a  previous 
exercise,  laying  doubt  about  the  trustworthiness  of  the 
w'ave  data.  The  w  ave  periods  and  wave  heights  were  not 
realistic  compared  to  buoy  observations  at  Sines  and 
nearshore  data  collected  near  Pinheiro  da  Cruz. 
Therefore,  we  do  not  show  any  comparisons  at  this 
location.  However,  we  feel  that  SWAN  tended  to  over¬ 
predict  wave  heights  in  this  general  area  possibly  due  to 
sheltering  effects  from  the  land  features  north  of  this 
area  for  which  we  did  not  account. 

6.4.  PC  Tides  versus  gauge  and  N6 

Comparison  of  PCTides  water  levels  from  tw'o 
different  data  sources  indicate  that  PCTides  provided 
very  good  agreement  in  amplitude  and  phase  during  the 
MREA04  Trial.  Fig.  10A  depicts  a  comparison  of  water 
level  from  PCTides  versus  data  collected  by  the  gauge 
data  collected  at  Sines.  PCTides  demonstrates  skill  in 
the  amplitude  and  phase;  the  correlation  coefficient  for 
this  comparison  is  1.0  w  ith  an  RMSE  of  9  cm. 

Fig.  I  OB  and  C  show-  a  comparison  of  w  ater  levels  at 
Pinheiro  da  Cruz  from  PCTides  versus  a  Nortek  current 
meter  deployed  during  MREA04.  The  mean  w  ater  depth 
from  the  period  of  record  from  the  Nortek6  (hence 
referred  to  as  N6,  sec  Fig.  5  for  location)  was  subtracted 
from  the  measured  water  depth  to  determine  the  water 
level.  Data  were  recorded  at  irregular  intervals  (approx¬ 
imately  every  60  or  120  min);  PCTides  water  levels 
were  interpolated  temporally  to  the  N6  observation 
times.  The  data  show  excellent  agreement  with  an 
RMSE  in  amplitude  of  8  cm,  phase  errors  less  than 
15  min  and  a  correlation  coefficient  just  under  1.0. 


Tabic  4 

Statistics  for  Deltl3D  versus  N6  for  period  April  4-11,  2004 


No.  of 
observations 

RMS 

Correlation 

Bias 

SI 

Hs  (m) 

103 

0.21 

0.70 

0.14 

0.33 

Direction  (°) 

103 

4.50 

0.70 

-  1.85 

0.02 

U-current  (m/s) 

103 

0.08 

0.53 

0.04 

-1.20 

V-current  ( m/s) 

103 

0.14 

0.66 

-0.04 

-1.35 

Water  level  (m) 

103 

0.07 

1.00 

0.02 
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6.5.  Delft  3D  versus  N6 

As  discussed  earlier,  Dclft3D  was  originally  run  in  a 
stepwise  fashion,  essentially  in  stationary  mode  for 
instances  in  time  at  three-hour  increments.  Later,  the 
configuration  was  adjusted  to  allow  for  a  ten-day 
hindeast  in  continuous  mode  and  included  an  improved 


configuration  and  an  updated  bathymetry  based  on 
newly  available  information.  By  the  end  of  the  exercise, 
data  from  all  the  in  situ  measurements  were  tabulated 
and  made  available  for  comparison.  Data  from  the 
center  gage  labeled  N6  were  compared  to  model  output 
of  currents,  wave  heights  and  wave  direction  at  the 
corresponding  grid  point  of  the  model  domain.  Analysis 


Fig.  11.  (A)  Time  series  of  DelfDD  model  output  and  N6  measurements  of  wave  direction.  (B)  Time  series  of  Delft3D  model  output  and  Nh 
measurements  of  significant  wave  height.  Tides  from  the  same  measurement  device  were  added  for  comparison.  (C)  Scatter  plot  comparing  model 
output  and  N6  measurements  of  mean  wave  direction.  (D)  Scatter  plot  comparing  model  output  and  N6  measurements  of  significant  wave  height 
(F)  Time  series  of  Delft3D  wave  height  with  and  without  winds  at  N6. 
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Measured  Wave  Direction  (degrees  nautical) 

Fig.  1 1  (i continued ). 


ofDclft3D  output  to  locations  N4  and  N9  indicated  little 
model  skill.  Uncertainties  in  the  nearshore  bathymetry  at 
these  particular  locations  are  attributed  to  this  discrep¬ 
ancy.  Table  4  presents  a  summary  of  statistics  for  wave 
heights,  wave  direction,  cross-shore  and  longshore 
currents  and  water  levels  from  the  model  versus  data 
collected  at  N6. 

Fig.  1 1 A  depicts  the  time  series  of  wave  direction  for 
both  Delft3D  and  N6,  whose  comparison  shows  good 
agreement.  With  the  beach  almost  facing  west,  the 
incoming  wave  angle  centered  around  270°  implies  that 
rip  currents  will  generally  occur.  It  is  expected  that  wave 
direction  would  refract  to  almost  perpendicular  to  the 
shore  by  the  time  wave  trains  reach  the  surf  zone.  Still, 
even  a  small  inclination  of  wave  direction  could  effect 
local  circulations  significantly. 

In  Fig.  1 1 B  the  time  series  of  significant  wave 
height  from  the  model  and  N6  is  plotted.  In  addition, 
the  water  levels  arc  plotted  to  help  point  out  that 
during  periods  of  low  tide  when  the  waves  are 
dramatically  affected  by  the  bar,  both  model  output 
and  N6  measurements  show  low  wave  heights  as 
would  be  expected.  The  bar  causes  the  waves  to  lose 
much  of  their  energy  due  to  breaking.  However,  during 
high  tide,  waves  propagate  over  the  bar  and  break 
closer  to  the  coast.  Fig.  1 1 D  depicts  a  scatter  plot 
comparing  the  modeled  significant  wave  height  to  the 
N6  gage  including  statistics,  which  show  a  consider¬ 
able  positive  bias.  Many  of  those  points  contributing 
to  this  bias  occurred  during  the  period  April  5-7,  in 


which  time  the  boundary  inputs  were  reaching  1  m,  so 
that  the  model  output  for  wave  height  at  high  tide  was 
substantially  higher  than  the  measurements.  The 
migrating  bar  certainly  contributed  to  the  discrepancy. 
This  discrepancy  may  be  also  due  to  the  fact  that  the 
boundary  conditions  result  from  a  running  string  of 
interdependent  models  rather  than  from  a  source 
ground-truth  data,  which  is  expected.  For  every  step 
taken  to  model  the  ocean  waves,  another  factor  of 
inaccuracy  is  introduced.  The  original  source  of  the 
error  may  have  been  in  the  regional  SWAN  model  due 
to  its  weakness  in  handling  the  sheltering  from  the 
energy  by  the  land  features  to  the  north. 

It  is  useful  to  examine  the  relative  importance  of 
wind  waves  against  swells  at  N6  buoy.  If  the  swell 
energy  is  dominant  for  the  whole  study  period,  then  the 
importance  of  using  a  good  local  wind  wave  model  is 
reduced.  The  comparison  of  Dclft3D  results  between 
with-  and  without-wind  inputs  is  presented  in  Fig.  1 1 E. 
Wind  waves  are  found  to  be  increasingly  important  after 
April  8,  when  the  input  swell  direction  veers  further  to 
the  north.  This  is  certainly  due  to  the  effect  of  headland 
blocking  (sec  the  land  features  in  Fig.  1 ).  The  increased 
blocking  is  consistent  with  the  wave  angle  plot  as  shown 
in  Fig.  8B.  The  inaccuracy  (i.e.,  wave  angle  and  direc¬ 
tional  spreading)  of  WAVEWATCH  III  spectra  could 
have  also  significantly  affected  the  Delft3D  results. 
There  is  no  measured  offshore  wave  data  for  verifica¬ 
tion.  However,  the  agreement  of  SWAN  output  at  the 
Sines  buoy  in  the  south  (as  shown  in  Fig.  8)  tends  to 
indicate  that  WAVEWATCH  Ill  provided  good  input 
spectra. 

Comparisons  of  the  components  of  currents  along  the 
shore  and  across  the  shore  were  plotted  separately  in  a 
time  series  and  showed  some  agreement  as  depicted  in 
Fig.  12A-B,  respectively.  At  high  tide,  very  little 
current  is  induced  due  to  wave  breaking  at  the  bar  and 
the  model  demonstrates  that  as  expected.  At  low  tide, 
the  model  output  indicates  there  should  be  an  induced 
current  in  both  the  cross-shore  and  longshore  directions. 
The  longshore  component  compared  well  with  a  corre¬ 
lation  of  0.66,  but  the  cross-shore  currents  compared 
less  favorably,  where  the  correlation  was  0.53.  The 
current  from  the  Delft3D  output  is  depth  averaged,  but 
the  measurements  are  at  a  particular  depth.  Since  long¬ 
shore  current  is  relatively  homogeneous  over  depth, 
depth  average  current  reasonably  represents  the  reality. 
On  the  other  hand,  the  cross-shore  current  consists  of 
significant  undertow.  This  makes  direct  comparisons 
somewhat  unrealistic  since  the  variability  of  the  current 
field  wnth  depth  is  substantial,  even  reversing  direction 
w'ith  depth. 
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Fig.  11  ( continued ). 


6.6.  Qualitative  evaluation  of  Delft  3D  verses  video 
imagery 

Comparisons  between  model  and  measurements  at  a 
single  point  give  only  a  limited  indication  of  the  model 
performance  and  observations.  For  a  more  complete, 
observable  picture  of  the  area  of  interest,  video  imagery 
was  taken  at  the  beach  site  for  which  snapshots  and 
time-lapse  averaged  imagery  were  generated.  This 
imagery  was  processed  into  geographically  rectified 
images  for  ease  of  comparison  to  data  also  geograph¬ 


ically  referenced.  Fig.  13  shows  a  time-averaged  image 
taken  on  about  April  7  at  10  GMT  with  DelfL3D  model 
output  points  showing  vectors  of  current  superimposed 
over  the  image.  On  the  image  the  lighter  shaded  regions 
associated  with  w'hitccaps  and  bubbles  are  indicative  of 
breaking  waves  and  from  this  one  can  infer  the  presence 
and  effect  of  the  bar.  At  intervals  along  the  bar  front, 
breaks  in  the  bar  can  be  inferred  where  less  wave 
breaking  occurred  and  more  cross-shore,  seaward  cur¬ 
rent  should  flow.  The  comparison  suggests  qualitatively 
some  agreement  between  the  Del03D  currents  and  the 
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Fig.  12.  (A)  lime  series  of  Dclft3D  model  output  and  N6  measurements  of  cross-shore  eurrent.  (B)  Time  scries  of  Dclft3D  model  output  and  N6 
measurements  of  longshore  current. 


inferred  locations  of  the  flow.  Since  the  bathymetry 
measurements  were  not  recent,  effects  of  the  bar  are 
expected  to  be  different. 

7.  Summary 

A  real-time  wave,  tide  and  nearshore  prediction 
system  can  provide  valuable  information  for  planning 
operations  by  describing  nearshore  characteristics  in 
littoral  areas  that  include  the  location  and  generation  of 
rip  currents  and  surf  conditions.  The  system  described 


in  this  paper  addresses  a  nested  modeling  approach 
with  a  scries  of  telescoping  models  acquiring  increas¬ 
ing  resolution  and  smaller  geographical  coverage  as 
one  approaches  the  desired  beach  loeation(s).  In  this 
study,  a  multi-nested  wave  model  is  used  to  provide 
boundary  conditions  to  the  nearshore  L)clft3D  model¬ 
ing  system.  Reasonable  bathymetry  was  available  for 
deeper  water  (greater  than  15  m  water  depth); 
however,  we  utilized  LIDAR  bathymetry  collected 
4  years  earlier  for  the  nearshore  modeling.  Uncertain¬ 
ties  in  the  nearshore  bathymetry  due  to  numerous 
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Fig.  1 3.  Time-lapsed  imagery  with  Dclfi3D  output  as  an  overlay.  Red  arrows  depict  the  magnitude  and  direction  of  the  current.  The  coordinate  system 
is  in  UTM. 


storm  events,  causing  the  movement  of  a  prominent  The  PCTides  tidal  model  performed  very  well  in  this 

offshore  bar  made  it  difficult  to  expect  predictive  skill  study  as  demonstrated  with  comparisons  to  measure- 
in  this  environment.  ments  with  a  shallow-water  instrument  and  tide  gaiige 
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data  at  Sines.  The  accurate  water  levels  were  provided  to 
the  Delft3D  modules  to  adjust  local  bathymetry  due  to 
tides. 

Reasonable  comparisons  were  obtained  when  com¬ 
paring  SWAN  wave  height,  period  and  direction  when 
forced  by  any  of  the  three  meteorological  models 
addressed  in  this  study  to  the  Sines  wave  buoy,  located 
about  45  km  south  southwest  of  the  beach  experiment 
area.  Improved  skill  is  expected  with  meteorological 
models  run  at  higher  resolutions  to  properly  resolve 
mcsoscale  features  including  land/sea  breezes  and 
orographic  effects.  These  results  indicate  that  the  atmo¬ 
spheric  models  are  capable  of  providing  sufficient  wind 
inputs  into  nearshore  wave  models. 

A  nearshore  Mcteo  buoy  provided  a  glimpse  of  the 
accuracy  of  the  meteorological  models  used  in  this 
study;  however,  uncertainties  from  the  wave  height 
sensor  make  it  challenging  to  adequately  assess  the 
skill  in  SWAN  at  this  location.  Future  efforts  should 
utilize  a  network  of  nearshore  buoys  for  a  means  of 
forcing  the  nearshore  modeling  system  (c.g.,  buoy 
directional  wave  spectra)  to  evaluate  its  predictive 
skill. 

The  results  from  Delft3D  were  mixed  and  illustrate 
the  need  for  more  accurate  input  information.  One 
factor  that  might  have  contributed  to  inaccuracies  of 
the  model  results  was  the  uncertainty  of  the  spectral 
boundary  conditions  from  WAVEWATCH  111  which 
could  not  be  precisely  evaluated  due  to  lack  of  ground 
truth  data  right  at  the  boundaries.  Another  factor  was 
that  DeIft3D  nearshore  waves  and  currents  are 
certainly  strongly  affected  by  the  migrating  sand 
bars  which  were  not  updated  due  to  lack  of  data. 
Though  the  results  gave  fair  indications  of  the  overall 
surf  conditions,  more  work  is  needed  to  improve  on 
the  operation  of  Dclft3D  to  provide  more  accurate 
results.  We  anticipate  future  improvements  towards 
the  collection  of  more  accurate,  up-to-date  nearshore 
depth  information. 
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